arXiv:hep-ph/9808250v2 6 May 1999 


The shock waves in relativistic superfluid 


G.V. Vlasov * 
February 1, 2008 


Abstract 

We consider the discontinuities in a two-constituent relativistic 
superfluid. In the acoustic limit they degenerate into the hrst and 
second sound which are independent up to the second-order linear 
approximation. Inclusion of the quadratic deviations relates to the 
small-amplitude shock. Particularly we consider a plane shock at low 
temperature when the phonon excitations contribute to the normal 
constituent. So we found the generalization of the temperature incre¬ 
ment and acoustic wave velocity in relativistic superfluid. The fourth 
sound speed is also calculated. 


1 INTRODUCTION 


The first approach to the relativistic superfluid mechanics proposed by Is¬ 
rael in and Dixon 0] concerns with perfect fluids. The method is useful 
for particular calculations, and as a general model can be applied to rel¬ 
ativistic superfluidity. For, strictly speaking, the coupled constituents are 
not perfect fluids: any coupling results to deviation from ideality. Or: the 
absence of coupling on microscopic level implies thermo-isolation of the con¬ 
stituents II. Nevertheless, the attempt of taking into account the deviation 
from perfect fluid is not senseless. The further development of Khalatnikov 
and Lebedev 0 includes the interaction between the superfluid and nor¬ 
mal constituent. This principle developed independently by Carter p was 
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gained in recent works |^, Althongh the main attention was directed to 
the general formalism, including the equations of motion, rather than applied 
problems. Carter and Langlois have recently derived the hrst and second 
sound speed in a superfluid with a phonon equation of state of the normal 
constituent. This inspires us to discuss the shock wave propagation in a two- 
constituent relativistic superfluid. While the relativistic shock waves hrst 
considered by Taub [^] have been under detailed discussion |ll|, |T4 

and the shock waves in superhuid helium were also investigated, the 
present task bears a qualitatively new feature: we work not in the frames of 
perfect huid hydrodynamics, for it is impossible to split the conserving par¬ 
ticle number current into a conserving ’’superhuid” and ’’normal” parts [§ 
as it is in the Newtonian limit |Q. 

We use the natural system of units {h = cnght = 1) and the metric 
corresponding to the Minkowsky space with a metric tensor diag {—h -l--l-). 


2 THE CONSERVATION LAWS 

The equations of relativistic superhuid mechanics are determined by the La- 
grangian L, whose inhnitesimal variation is given by formula |^, |^, 

5L = (1) 

where the particle number vector conjugated to the momentum co-vector 
/ig obeys the conservation law 


Vgn^ = 0 (2) 

and the entropy vector conjugated to the thermal momentum 1-form 0^ 
is also conserved 

V,s^ = 0 (3) 

till a shock wave appears. The energy-momentum tensor corresponding to 
the Lagrangian L has the form 

(4) 

with the pressure function 

T = L - (5) 
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( 6 ) 


The conservation law of the energy-momentum tensor 

V,T,^ = 0 

leads to the equation of motion of the normal constituent 

V.0,) = 0 

and the irrotationality condition 

^ g ^ u 0 


(7) 


( 8 ) 


Vectors and /r^(also and 0^) are not colinear. For the particle number 
vector determines the Eckart rest frame, while determines similarly 
the superfluid rest frame and that does not coincide with the former. The 
Lagrangian L (/i, y, s) depends on three invariants, them being 


§2 = -S^Sr. 


(9) 


for the normal rest frame entropy density; for the cross product given by 

1/2 = ( 10 ) 

and for the effective mass variable (i.e. the chemical potential in the super¬ 
fluid rest frame) 

( 11 ) 

The secondary variables and 0^ can be expressed through the primary 
variables and according to the formula 




B -A 
A C 




( 12 ) 


with the coefficients obtained immediately by differentiation of the La¬ 
grangian 

On the other hand we can write 




F Q 
Q G 


QQ 


B + A^jC -AjC \ ( 
-A/C 1/C ] i 
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(14) 




and calculate new coefficients through the pressure function 


PI 




( 15 ) 


However, a more convenient way is to obtain all parameters ([^) in terms 
of (D, diD and (|TI|) and express them in terms of the relative translation 
speed between the superfluid and normal reference frames 

= 1-— (16) 

and the effective temperature 


02 = (17) 

instead of invariants y and s. At low temperature the phonon-like excitations 
with energetic spectrum uj = cp, where 

2 ndjji 
/i dn 

and the latter is the sound speed contribute to thermodynamical functions 
of the normal constituent. The ’’phonon” Lagrangian has the form [§] 

2/3 ^ 0 " 

90c3 [1 — 

where the pressure of excitations is 



L = P — xjj = 


n/3 


s2 + 


0 - 1 ) 





and 


4 

g _ 0^ 

\/l — w'^ 15c3 [1 — t(;2/c2]^ 


( 20 ) 

( 21 ) 


is the entropy g^ in the superfluid reference frame, while fc is a dehnite con¬ 
stant. The coefficients (P) then are 


1 -c2 0 

/i c2 — t(;2 


c 


1 e 

Ss c2 — r(;2 


( 22 ) 


For an ultra-relativistic spectrum of excitations (c —> 1) the term A disap¬ 
pears. 


4 



















3 THE DISCONTINUITIES 


The discontinuities in superfluid helium have been first discovered by Kha- 
latnikov [§. Discussing the discontinuities in a relativistic superfluid we shall 
follow the standard formalism of the relativistic shock waves [^. However, 
since we deal with a strongly self-interacting medium and two-constituent at 
that, a standard perfect fluid theory is impossible. So we have to think 
of a new extended formalism, combining the theory of relativistic shock 
waves 0,0 P!3| , and the relativistic superfluid mechanics |§, 0, § • In the 
Newtonian limit this method must be reduced to the non-relativistic theory 
of discontinuities in helium |^, and in the acoustic limit this method must 
give the first and the second sound in a relativistic superfluid PJ. The am¬ 
plitude of discontinuities is assumed to be not very large, since superfluidity 
is expected to take place on both sides of the front with no phase transition. 

Let the hypersurface S be the front of the discontinuity and vector 
be a unit (space-like) normal to it. The conservation laws (2), (|^) entail the 
conditions 

|A,n«] = 0 [X,T;] = 0 


i.e. 

= A^n^ (23) 

\,Tl, = A,Tf, (24) 

where indexes -|- and - relate to the quantities ahead of and behind the 
discontinuity, respectively, and the square brackets imply the change across 
the front of a discontinuity. Eq. (p3|) conveys the continuity of the orthogonal 
part (marked always by index T) of the particle number current. 

We cannot make use of the entropy conservation law (^) since it does not 
take place in the shock waves. For that reason the results of Ref.[0 pertain 
for a most general two-constituent system (but not a superfluid itself), while 
the discussion of small perturbations suffers no restrictions concerning eq. (|^). 
Indeed, the results for shock waves in a two-constituent system and in a 
superfluid will coincide when the magnitude of the relevant discontinuities 
tends to zero. 

For a two-constituent superfluid the additional equation is the irrotation- 
ality condition (|[) which yields 
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Multiplying it by and, then, by a unit vector orthogonal to A^, we get 
an important relation 

A = A (25) 

(where index || denotes the tangential part of any quantity, particularly /i^ = 
This is the very condition for conservation of irrotational motion 
passing a plane shock wave [TO. It is clear that for a multi-constituent 


system the irrotational motion conserves at both sides of the shock wave if 
condition (j^) takes place for each constituent. 

Substituting the expression (^) and (^) in we obtain 


n~ 


- iJ-„) + sie+„ - =-(*+- S-) A, 


(26) 


Multiplying this equation by A'^, we get 

A+ ~ A) + -5+0+ - = - (\k+ - T_) (27) 

Then, multiplying eq.(^) by the unit vector rj’^, in light of (P5|), we get 

(^A ~ A) + s:[:0+ — 51(01 = 0 (28) 

This important relation determines the types of discontinuities. Then, the 
irrotationality condition ( [25|) than coincides with the condition of strong 
discontinuity in the superfluid constituent which occurs when 7 ^ 0 . As for 
the discontinuity in the normal component, in light of (^^, it is determined 
merely by the single relation 


5701 = S70i 


(29) 


Thus the constraint (^) corresponds to ordinary shock waves, while (p^) 
beseems to a ’’temperature” discontinuity of the second sound type. 

As a particular instance of another type of discontinuity we consider a 
vortex sheet in superfluid |]^. Since the vortex sheet separates the whole 


space into domains where the superfluidity takes place, the irrotationality 
condition (|]) does not hold in the global sense and we cannot establish the 
constraint (P5[ ) at both sides of the sheet. Therefore, the tangential discon¬ 
tinuities are possible. The conservation law (^) then yields = n(: = 0 
that determining weak or slip-stream discontinuity: the particle number flow 
across the front of discontinuity equals zero, indeed, no matter crosses the 
hypersurface of the discontinuity, i.e. this hypersurface is made up of stream¬ 
lines of the fluid. 
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4 A PLANE SHOCK IN FLAT SPACE 


Let the discontinuity propagates along the axis . We choose the unit normal 
= (0,1, 0, 0) - and the medium at rest before the front. As a rule one used 
to practice with the rest frame co-moving the front, so that the fluid flows in 
the front with the velocity which is equal to that of a shock wave. Hereby, 
the relevant vectors and co-vectors may be presented as 

n^ = n(^^l + (p2,(p,0,o) pg =/i (^—yiT^,^,0,0^ (30) 

= s(Vl + a2,a,0,0) 0, = 0 A 0, o) (31) 

Since the medium ahead of front is at rest, the relative velocity equals 
zero, = /i-S-, while behind the shock 


yl 




VT 




Also 


a_ = /3_ = (p_ = ^_ = a: 

Hence, the velocity of the shock is determined as 

X 


u = 


ViT 




(32) 


(33) 


Substituting our dehnitions 
we, hrstly obtain 


and (^) in eqs. (|^), (^5]) , (|^) and 


n+(p = n_x (34) 

The rest equations, in view of ([5^), will be 

=/i-VTT^ (35) 

05+0+^! + (3"^ = xs_0_a/1 -|- x"^ (36) 

ri-X (^/i+ — x/i_) -|- a/3s+0+ — x^S-0- = —4/+ -|- (37) 


The parameter (p incorporates only in eq. (^) and it can be calculated as 
soon as the rest unknowns are found. Thus, in eqs. (^), (|36|) and ( 0 ) 
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the unknowns are: the four parameters a, f3, x and three invariants (j^), 
(p!0D , ([TT|) behind the shock on which the pressure \I/+ depends. The pressure 
behind the shock can be expressed through /i+, 0+, and through the relative 
velocity for which we use the notation w. The formula 


= Vl + 1 + 


(38) 


Vi — 

relates the later quantity with a and Our goal is to hnd the velocity of 
the shock wave u for a single parameter given behind the shock. Without 
the loss of generality can be chosen for this parameter. Thus, there are 
six unknowns in four equations (^51) , (^6]), (^), (^). The rest two relations 
follow from (|T^) or ([T^ ) 


n_x = F+/i+^ + (5+0+/5 


(39) 


s_|_q; — T G+0+/3 


(40) 


with the coefficients (p!3D calculated for the state behind the front. 

The knowledge of the equation of state in explicit form is necessary for 
calculation of the right-hand side of eq. ( 0 ) and the coefficients in eqs. (p9|) 
and (M). 


5 A LOW TEMPERATURE CASE 

The low-temperature equation of state was derived by Carter and Langlois 


|]. In view of the expressions 

form 

'i/’x/l — w‘^a\jl + = X Vl + x"^ 


, (^0[) take the 
(41) 


X — x) + T — w'^ — 4x'^ + ‘ip — ^ (^P — 


/iF^ + (l - 4) r^/3 


= X 




\/l — w'^a = (l — 4) ^ + (|4 ~ 


(42) 

(43) 

(44) 
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where 

i’- 


fi-n_ 


r = 


/i_n_ 

P_ 


F 



At low temperature we have the estimations 


F = 


n 


/i 


0(0^) 


G ~ 02 Q ~ 0^ 


( 45 ) 


implying that t = O (0^) and, hence, equations (^Tf ) and (|3^ approximately 
(up to the terms O (0^)) coincide with their zero-temperature version. 


6 THE SOUND, STRONG AND SMALL- 
AMPLITUDE SHOCK WAVES 


If all parameters behind the front tend to their values ahead, the shock be¬ 
comes a sound wave. Since the entropy in the sound wave is conserved, we 
can apply formalism achieved for a two-constituent relativistic medium 
with the conserved particle currents of both constituents. In the linear ap¬ 
proximation both methods lead to the same result, namely from the system 
(|35D, (^), (|37|) , (pSD, (p^, (|40|) we obtain an equation for two branches of 
sound at arbitrary temperature which is analogous to that derived by Carter 
[|] and, under assumption = o{C), splits into 


u^i = 


-B 


B -|- fiB^ F sA^ 


(46) 



sBg + fiAg 

^ C 


(47) 


and reduces, in the low temperature limit, to the hrst and the second sound 
speed, respectively ^ : uj = c, un = c \/3. Here for any variable V we used 
the notation 


ds 2y dy ^ dfi 2y dy 


(48) 


However, if we omit A, the second sound speed calculated by formula (^) 
with the phonon Lagrangian (^) of Carter and Langlois p[ will be uu = 1/3 
instead of obvious uu = c/3. Because the Lagrangian (|^) is derived for the 
two-fluid theory with non-zero cross term; while the Lagrangian of thermal 
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excitations of the Israel theory |]^ differs from ([I9|); although the relative 
translation speed w between the constituents is presented in both approaches. 
So each Lagrangian is useful in the theory to which it does belong. 

In order to hnd the velocity of a small-amplitude shock we rewrite 
eqs.dD, (HI), (HD, (P, (H), (H in the second-order approximation. After 
tedious calculations we hnd the velocity increment 


An 


II 


A0 


Uii 


1 + C2 0 


(49) 


of the shock corresponding to the second sound when the temperature in¬ 
crement A0 tends to zero. In turn, the latter gives rise to a hnite relative 
speed w behind the shock, since 


A0 


w 




(60) 


In the non-relativistic limit the equations (HD (H) yield well known ex¬ 
pressions 0. The shock occurs ahead of the second sound, as is in superhuid 
helium at low temperature. 

So, in the acoustic limit the solution splits into two independent branches 
Ui and uji corresponding to the hrst and the second sound. The hrst branch 
describes wave propagation through the medium which behaves as a perfect 
fluid composed of two constituents whose pressure and enthalpy are P -|- "0 
and /in -|- Aip respectively. While the constituents in the waves of the second 
branch move independently, a counterflow appears: tc 7 ^ 0. In general, a 
’’mixed” solution occurs, and the temperature increases together with the 
chemical potential. 

The estimations (HP imply that the hrst sound and the relevant shock 
wave coincide roughly with usual discontinuity in the cold constituent.A great 
pressure jump AT is produced inevitably by the change in pressure of the 
superhuid constituent P, since the contribution of the normal constituent 'ip 
is small. Hence, in view of eqs. (H) and (H) we conclude that strong shock 
waves at low temperature propagate with the speed u = uq + O (0^) which 
approximately equals to the speed of a usual shock wave in cold constituent 
Mo, hut always ui > uu- A more precise result is 


Ur = u, 


{l-rul (4-r)} 


(51) 
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Table 1. The gamma-factor 7 ^ = w/^Jl — w"^ of the relative translation speed 
vs the Mach number M = ujcoi the shock wave and the pressure change 


1 

1 

+ 

1.1 

1.5 

2 

5 

10 

M 

1.024 

1.10 

1.17 

1.35 

1.46 

'^W 

0 

1 

CO 

.031 

.065 

.256 

.476 


^ 7 / = Mg {1 “ ™o (4 — r a/I — I (52) 

The second-sound discontinuities should be regarded as ’’moderate” for 
intermediate values of «//. This takes place if relative changes in the super¬ 
fluid and normal variables are of the same order and they can be of the same 
order if they do not access r greatly. For a superfluid matter of neutron 
stars and phonon equation of state it is easy to estimate r 10 For this 
particular example we performed calculations with an ultra-relativistic su¬ 
perfluid matter. The sound speed in this medium equals exactly to c = 1/a/3, 
and approximately it is the hrst sound, while the speed of the second sound 
cii = c/ a/3- For a not very small pressure change the velocity of the shock 
wave will be merely Uq. The dependence of the relative translation speed w 
on uq is given in table 1 . 

Although the second sound velocity attains to the saturation value mq, 
the relative translation speed w grows with the growth of the shock wave 
intensity. 


7 THE FOURTH SOUND 


The fourth sound takes place when the normal constituent is restrained by 
some external agent, while the sound propagates through the superfluid con¬ 
stituent. We cannot use the equation of motion (|^, but the conservation 
laws (H), (H) and the irrotationality condition (||) will determine the fourth 
sound speed. If the sound wave propagates in the direction determined by 
vector = (u, 1, 0, 0) the change of gradient of arbitrary quantity V is pro¬ 


portional to its inhnitesimal change V : 


VW 


= VXv. Thereby, we write 


— uh^ - 1 - = 0 

-ujxi - /io = 0 
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(53) 

(54) 














- MsO + = 0 (55) 

instead of (^, (^, and (|^. In the reference frame co-moving with the 
normal constitnent we pnt 


= s (1,0,0,0) /ig =0,0) (56) 


where /i„ = /i/ a/I — w"^ is the chemical potential in the ’’normal” refer¬ 
ence frame. Since the discontinnities propagate throngh the snperflnid con¬ 
stituent, there must be 

= 0 (57) 


Equations (b^, (^J), (|^ analogous to the relevant non-relativistic set [Q 
determine the speed of the fourth sound. Requiring the vanishing determi¬ 


nant of the system (^^, (|5^) we get a quadratic equation for u. While the 
speed of the hrst and the second sound is determined by a 4-order system [|], 
the fourth sound speed follows from two equations. At zero temperature the 
speed of the first and the fourth sound are obtained by the same equations 
(1^ and (^) and coincide exactly with the sound speed c in the superfluid 
constituent. The difference appears at hnite temperature on account of the 
relationship g] 


between the infinitesimal discontinuities in ([53|), (^) and the temperature 
dependence of matrices in (|^). For the phonon equation of state we get the 
explicit formula 


M/y — C -j- C 


Pn 

Ps 



-3c^ - 


cp 


dc 


dh 




which generalizes the relevant non-relativistic relation |^, where pn and ps 
is the normal and the superfluid energy-density respectively |@|. 


8 CONCLUSION 

Summarizing the results obtained in the present study, we emphasize the 
formulae (0), 0), 0, (0 which determine the propagation of disconti¬ 
nuities through a two-constituent relativistic superfluid in the general case. 
For a plane shock wave seven equations (^), (^), (|S6|), (|T7D , (|3^, (|3U[) and 
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(^0[) with seven unknowns must be solved. In the acoustic limit these equa¬ 
tions reduce to formulae (^), (^71) for the hrst and the second sound. At low 
temperature the system (pj]), (|55D, (PS|), (P7|), (|55D, ( PP|) and (^01 ) reduces 
to (|34D , (^, (^ , (^ , (^, and (^. The velocity of strong shock 
waves is given by (^) and (^), while and describ the change of 


parameters in a weak shock wave. As for perspectives and applications, the 
shock waves and spin-isospin sound in the nuclear matter are worth to be 
discussed in future. 
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